Absent pinch points and emergent clusters: further neighbour interactions in the 

pyrochlore Heisenberg antiferromagnet 
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We discuss the origin of spin correlations observed in neutron scattering experiments on the para- 
magnetic phase of a number of frustrated spinel compounds, most notably ZnC^CU. These corre- 
lations are striking for two reasons. First, they have been interpreted as evidence for the formation 
of weakly interacting hexagonal clusters of spins. Second, they are very different from those cal- 
culated for the nearest neighbour Heisenberg pyrochlore antiferromagnet, in which Coulomb phase 
correlations generate sharp scattering features known as pinch points. Using large-n calculations 
and Monte Carlo simulations, we show that very weak further neighbour exchange interactions can 
account for both the apparent formation of clusters and the suppression of pinch points. 

PACS numbers: 75.10.Hk, 75.50.Ee, 61.05.fd 



I. INTRODUCTION 

A central reason for interest in geometrically frustrated 
magnets is that they are systems which can support 
strong correlations without developing long range or- 
der. This cooperative paramagnetic behaviour appears 
at temperatures small compared to the main interac- 
tion scale, which is characterised by the magnitude of 
the Curie- Weiss constant 0cw- The absence of long 
range order at low temperature can be understood in the 
framework of classical models with nearest neighbour in- 
teractions, as a consequence of macroscopic ground state 
degeneracy— 

Results from recent neutron scattering studies of some 
important examples of geometrically frustrated magnets 
present a paradox, and our aim in this paper is to offer 
a resolution of the paradox. More specifically, experi- 
ments on a number of spinel compounds, including in 
particular ZnCr204, have been interpreted as revealing 
the formation of small, independent clusters of spins*—. 
In contrast, the minimal theoretical model which is as- 
sumed to describe these systems, the nearest neighbour 
Heisenberg pyrochlore antiferromagnet, is notable for a 
quite different feature: pinch point scattering, indicative 
of algebraic correlations. We show in this paper that fur- 
ther neighbour interactions strongly influence scattering 
in the cooperative paramagnetic phase. Very weak fur- 
ther neighbour interactions of the appropriate sign (fer- 
romagnetic second neighbour, or antiferromagnetic third 
neighbour) are sufficient to suppress pinch-point scatter- 
ing and generate quite accurately the short-range correla- 
tions that have been interpreted in terms of independent 
spin clusters. 

The spinel compounds, AB2O4, have a structure in 
which the S-sites form a three-dimensional network of 
corner-sharing tetrahedra, known as a pyrochlore lattice 
(Fig. [1]). The same lattice is realized by both A and B 
sites in the pyrochlore compounds A1B2O7, and in the 
Laves phase compounds with formula AB2 . Evidence for 
spin clusters has been reported in the chromites, ACr 2 04 
with non-magnetic A = Zn, Mg and Cd as well as in the 



ferrite CdFe2C"4 (Refs. [2|-|5| respectively). In all cases, 
the strong similarity between experimentally measured 
spin correlations and the form factor for independent 
hexagonal clusters leads to an interpretation in terms of 
independent spin clusters. Results for Y(Sc)Mn2 have 
broadly similar features-. Cluster like correlations are 
also observed in the spin ice compound Dy2Ti20y: here 
it has been argued from comparison with simulations that 
they arise from further neighbour interactions^. Spin 
clusters are not universally observed in frustrated py- 
rochlore magnets, however. Among the spin ice com- 
pounds, compelling evidence for pinch point scattering is 
seen in Ho2Ti20v^. For a recent review of properties of 
the frustrated spinels, see Ref. |3 

The classical nearest neighbour antiferromagnetic 
Heisenberg model on this lattice (which we refer to as 
the NN model) is now well understood theoretically-^— . 
The model has no phase transition and remains in the 
paramagnetic phase down to temperature T = 0. The 
low temperature regime is nevertheless strongly corre- 
lated. The nearest neighbour interaction energy is min- 
imised by any spin configuration for which the magneti- 
zation on every tetrahedron vanishes, giving a ground 
state manifold with a macroscopic number of degrees 
of freedom. While states on this manifold are weighted 
differently by thermal fluctuations, order-by-disorder ef- 
fects are too weak to induce a Neel state. Simulations 
of neutron scattering show sharp features in the struc- 
ture factor—. Known as pinch points, a clear under- 
standing of their origin comes from the recognition that 
ground states of the NN model can be mapped onto con- 
figurations of a solenoidal lattice vector fiel d 14 ! 15 (this 
mapping will be reviewed in Sec. IIIBI) whose long wave- 
length fluctuations are governed by an entropic Maxwell 
action. The pinch points in scattering are a consequence 
of dipolar correlations of the solenoidal field. (For a re- 
cent review of models with such a Coulomb phase, see 
Ref. [lfl) Although algebraic at T = 0, at finite tem- 
perature spin correlations decay exponentially beyond a 
correlation length £11, and pinch points acquire a width 
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FIG. 1. (color online) Pyrochlore lattice, illustrating the first, 
second and third neighbour couplings. There are two inequiv- 
alent third neighbour interactions, but for J^ a ~ Jsb, we use 
J3 to denote both. Numbers indicate sublattice labeling. The 
centres of the tetrahedra form a diamond lattice and the red 
arrows show the choice made to orient the diamond lattice 
bonds. 

~ f 7 . For spin ice, the correlation length diverges ex- 
ponentially as T —> due to the finite energy cost for 
violations of the ice rules, but for the NN model with 
continuous degrees of freedom, only as T" 1 / 2 . Since the 
NN model unambiguously gives rise to pinch points in 
scattering when T <§; 0c w, with thermal broadening at 
finite temperature, their absence in experiments on the 
paramagnetic phase of the frustrated spinels demands an 
explanation which goes beyond the NN model. 

Here we extend the NN model by including further 
neighbour exchange terms (see Fig.[T]for definitions) and 
consider their effect on paramagnetic scattering. Fur- 
ther neighbour terms, both ferromagnetic (FM) and an- 
tiferromagnetic (AFM), will in general have two effects: 
(i) to change scattering in the paramagnetic phase; (ii) 
to induce ordering at a low T. With regard to or- 
dering transitions out of the low temperature param- 
agnetic phase, the role of further neighbour exchange 
has been considered before: from a mean field perspec- 
tive, Reimers et al— examined possible ordered states 
with further neighbour exchange; and following simula- 
tions by Tsuneishi et al. 18 , a detailed study by Chern 
et al— of the phase diagram with a weak J2 confirmed 
the mean field predictions. There has been less work 
to pin down the effect of further neighbour exchange on 
paramagnetic scattering in Heisenberg models, although 
we note that an analysis of neutron scattering data us- 
ing third neighbours has been done for the frustrated 
spinel ZnFe2 04, which has ferromagnetic nearest neigh- 
bour interactions 20 . In spin ice, Yavors'kii et al— have 



shown with extensive simulations that further neighbour 
interactions can promote cluster-like scattering. 

Our main results are as follows: first, further neighbour 
exchange terms of the right sign (FM J2 or AFM J3) can 
be effective both at suppressing pinch point scattering 
intensity and at broadening pinch points; and second, 
within the manifold of NN ground states, they favour 
states with cluster-like correlations. These effects are 
noticeable even at temperatures two orders of magnitude 
larger than the temperature of the transition that the 
additional interactions induce. Using the language of the 
Coulomb phase description, these effects are summarised 
in part by the fact that further neighbour terms renor- 
malise the Coulomb phase coupling constant. 

A visual impression of our findings in relation to ex- 
perimental reports of spin clusters is given by compar- 
ing figures that show the computed diffuse scattering in 
two cases. The first case is the nearest neighbour model, 
for which pinch points, broadened by finite temperature, 
are apparent in Figs. [3] (a) and [3] (b). The second case 
is a model with additional weak further neighbour in- 
teractions, for which behaviour is shown in Figs. [3] (c) 
and UJ (d). In this case scattering at the positions of 
pinch points is heavily suppressed. In its place, rings 
of scattering intensity are formed of the kind that have 
been interpreted as indications of spin cluster formation 
in Refs. SHl 

The remainder of this paper is organised as follows. In 
Section |TI] we present the model and detail the main ana- 
lytic tool used, a self-consistent Gaussian approximation 
which is exact for n-component spins in the large-n limit. 
Here we also review the mapping to lattice vector fields, 
which facilitates a more straightforward discussion of the 
long wavelength physics in the vicinity of pinch points. 
In Sec. IIIII we consider the effect of a FM J2 or an AFM 
J3 which are of the sign required to induce cluster-like 
scattering and suppress pinch points. In Sec. IIV1 we 
summarise the effects of AFM J2 or FM J3 , and we look 
briefly at other combinations in Sec. [Vj In Sec. I VII we 
discuss dynamics, describing how we expect the effects 
of further neighbour interactions to be more conspicuous 
in quasi-elastic scattering than in static correlations. We 
consider estimates of further neighbour exchange for rel- 
evant compounds in Sec. IVIII We summarise our results 
in Sec. I Villi Further details of various calculations can 
be found in Ref. [H 

II. MODEL AND THEORETICAL APPROACH 

Our focus is on the low temperature behaviour in the 
paramagnetic phase of the classical isotropic Heisenberg 
model with Hamiltonian 

// i^\/,„S,.S ( (I) 

where the degrees of freedom are unit length, three- 
component vectors residing on the sites of the pyrochlore 
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lattice and where the dominant interaction is nearest 
neighbour antiferromagnetic (see Fig. [IJ. It is known 
that behaviour of the nearest neighbour Heisenberg 
model is very well approximated by the large-n limit of 
its generalisation to n-component spin a i 22 . 

We also use this large- n limit. Through its equiva- 
lence to a spherical model 23 ', we present our method in 
the language of a self-consistent Gaussian approximation 
(SCGA). The validity of this approach is confirmed by 
the excellent agreement with Monte Carlo (MC) simula- 
tions on the Heisenberg model, Eq. (TTJ), both with and 
without further neighbour terms - results from Monte 
Carlo simulations appear as points in figures and results 
from the following analysis as solid lines. Our simula- 
tions were performed using the Metropolis algorithm for 
systems of 2048 spins, and 16384 spins. The longest runs, 
required for the lowest temperatures, were 10 8 MC steps 
per spin. In general though, run lengths were chosen 
in order that the statistical errors are smaller than the 
symbol size in the plots. 



A. Self-consistent Gaussian approximation 

We replace fixed length spins Si with soft ones, de- 
noting the magnitude of one component by Si . Soft spin 
configurations are weighted by e -/3w where 



m = \ E( A % + /?E Jnvtfrwj 



(2) 



The value of the Lagrange multiplier A is chosen to en- 
sure (sf) = 1/3 as required for a single component of 
a unit length three-component vector in a paramagnetic 
phase. The m-th neighbour interaction matrices Vfj 
are related to the respective adjacency matrices by ad- 
dition of multiples of the identity in the way described 
in Appendix [B] and these (seemingly arbitrary) choices 
of the diagonal entries are dictated by our subsequent 
interpretation of A as the stiffness of a Coulomb-phase 
action. The multiplier A is determined self-consistently 
by the condition on the average spin length 
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where AN is the total number of sites. The translational 
symmetry of the system ensures that the interactions are 
block-diagonal in reciprocal space (see Appendix [A] for 
Fourier transform conventions), and the self-consistency 
condition can be written in that basis as 
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AI+/3^J„VW(q) 



(3) 



where the matrices \A")(q) act within the space of the 
four sublattices. Their explicit forms are provided in 



J 3 = O.OUi 
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FIG. 2. Dependence of flux stiffness A on temperature, show- 
ing influence of weak further neighbour interactions. Solid 
lines: solutions of Eq. (O for the cases J3 = ±0.01Ji and 
J3 = 0, all with J2 = 0. Points: results from simulations on 
2048 spins interpreted via Eq. @. The effects of third neigh- 
bour interactions are clearly noticeable even at temperatures 
two orders of magnitude larger than J3. Insets: detailed be- 
haviour at low temperature; axes as in main panel. First order 
phase transitions occur in simulations for both signs of J3. 



Appendix [B] The solutions for A for various choices of 
further neighbour exchange are shown in Fig. [2] We dis- 
cuss analytical bounds on the solution following Eq. (|f 01) , 
after the physical significance of A as the stiffness of a 
Coulomb-phase action has been made clear. 

The correlation functions in reciprocal space between 
different sublattices arc 



Mq)4(q)> = 



AI + /3^J„V(")(q) 



(4) 



J \1V 



and the static structure factor is the sum of all 16 such 
correlators. For the nearest neighbour model this can be 
obtained analytically at all temperatures in a palatable 
form: at T = 0, it is provided by Isakov et al.— (with 
a different sublattice labeling convention) and for com- 
pleteness we provide the full, finite temperature form in 
Eq. (IDip in Appendix [D] In the general case where fur- 
ther neighbour terms are not all zero and analytic expres- 
sions are unilluminating, spin correlations at any temper- 
ature can be straightforwardly obtained numerically from 
Eq. Q. 



B. Mapping to long wavelength description 

Within the SCGA, Eqns. @ and © give a full de- 
scription of spin correlations in the paramagnetic phase. 
A more physical understanding comes from studying the 
long wavelength physics in a way that allows comparisons 
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with related models, including spin ice, interacting dimer 
models, and more generally any system with a Coulomb 
phase, such as those discussed in Ref. [TJ| and references 
therein. 

In our context, the long wavelength degrees of free- 
dom emerge naturally after interpreting spin configura- 
tions as configurations of a lattice vector field defined on 
the bonds of the diamond lattice. The mapping—^ is 
constructed as follows. Each site on a pyrochlore lattice 
is associated with a diamond bond, whose midpoint it is; 
the diamond lattice is bipartite so its bonds are orientable 
in a consistent way across the lattice. Therefore each spin 
component can be interpreted as the local magnitude of 
a lattice vector field directed along diamond bonds. The 
merit of this description is that the vector field obtained 
is divergenceless in the limit /3J± — > oo at fixed j3J n >i- 

More precisely, for each spin component, labeled by a, 
define the flux field B Q by B Q (r) = S a (r)e(r) where e(r) 
is a (real space) vector oriented along the diamond bond 
passing through the pyrochlore lattice site r. The map- 
ping has a straightforward implementation as a change 
of basis in reciprocal space: 



/M Q (q)N 
5, Q (q) 
^ Q (q) 




(5) 



The directions of the vectors e(r) and the sublattice la- 
beling are shown in Fig. [1] The vectors e(r) are taken 
to have length \/3/4 (in contrast to the convention of 
unit length by Isakov et al.—) so that the definition of 
the mapping is consistent with the orthogonal change of 
basis in Eq. ([5]). At q = 0, the three spatial components 
of B are the three staggered magnetization order param- 
eters which would describe a transition to q = Neel 
order as referred to in Refs. QjJ and Hi]. The experimen- 
tally accessible quantity is the structure factor which, we 
define as S(q) = 4(|M(q)| 2 ). 

We are interested in the small q behaviour of the vec- 
tor field B(q) but the experimental probe is via the linear 
combination M(q). Consequently, the long wavelength 
physics does not appear near q = in the structure fac- 
tor, but rather near q = K for special vectors K, which 
satisfy M(q + K) cx -Bj(q) for some spatial direction i. 
The K are a subset of the reciprocal lattice vectors (see 
Appendix IC|) . 

Using this mapping, we now turn to the long wave- 
length physics of the model defined in Sec. Ill Al We con- 
sider explicitly the case of J^ a — J 3 f, = J3 with all other 
further neighbour interactions zero; it is straightforward 
to carry out a similar analysis with any other choice of 
further neighbour terms. 

At small q and j3J\ S> 1, the linear combination M(q) 
is nearly frozen out and within the SCGA we can le- 
gitimately read off a small q theory for the field B by 
restricting the interaction to the 3x3 subblock which 
acts between the components of B, providing A is still 



calculated correctly using the full theory, Eq. ([3]). Fur- 
ther details are provided in Appendix |Bj By expanding 
the interactions to second order in q we obtain as a long 
wavelength theory 

m = \T. A l B (q)l 2 + P J i a2 \i ■ B (q)l 2 

q 

-8/3J 3 aV|B(q)| 2 . (6) 

The length scale a is the nearest neighbour distance on 
the pyrochlore lattice. It is clear from Eq. © that J\ 
acts purely to enforce the constraint q • B = 0. 

We now consider the flux correlations that follow from 
Eq. ©. The q = flux fluctuations are straightforwardly 
obtained as 



1 



(A(q)£*(q)) q=0 = -<5. 



(7) 



from which we identify the parameter A as the dimen- 
sionless stiffness of the flux fields in the long wavelength 
description. Since its original role as a Lagrange multi- 
plier is to imitate the restrictions on phase space due to 
the fixed spin length, it is natural that it appears also as 
the stiffness which is entropic in origin. The correlation 
function in Eq. © is directly accessible from simulations 
and the data points in Fig. [2] follow from this relation; 
the same relation is used to measure the stiffness in sim- 
ulations of dimer models (see e.g. Eq. (2) in Ref. H^). It 
is impressive that the SCGA captures behaviour in the 
Heisenberg model so accurately. 

The correlator for longitudinal flux fluctuations is 



(B x (q x )B*(q x )) qy=qz=0 = - 



(8) 



where we introduce the longitudinal correlation length 
£|| = a \J /3(Ji — 8Jz)/\ (whose explicit form of course 
depends on which further neighbour exchanges are in- 
cluded). The inverse correlation length is the width of 
pinch points in reciprocal space. In the limit — > 00 
at fixed /3J3, flux fields are strictly solenoidal and £m di- 
verges. More generally, we emphasise that £y depends on 
P both explicitly and through the /3-dependence of A, via 
Eq. ©. 

The transverse flux correlations are 



(B z (q x )B* z (q x )) 



I 



-, _ 8/jJ 3 a 2 „2 
1 \ 



(9) 



For J 3 negative, the length £1 = ay/8(3\ J3I/A is the 
transverse correlation length for flux fluctuations. For 
positive J3, £± is not a correlation length: instead it sets 
a length scale at which higher order terms in q must be 
included in the long wavelength theory. 

Youngblood and Axe— considered a phenomenological 
Landau free energy for polarization fluctuations in ferro- 
electrics, essentially identical to Eq. ©. Our discussion 
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differs in that Eq. ^ is derived directly from a micro- 
scopic model, the parameter A is allowed to flow, and 
moreover both positive and negative coupling to trans- 
verse gradient terms are considered. 

Having established that A has a physical interpretation 
as the Coulomb phase coupling constant, we now turn to 
bounds on the value of A which solves Eq. (|3|) . Let e m i n be 
the minimum eig envalueof £ n J„VW(q) over all q. The 
requirement that the quadratic form in Eq. @ is positive- 
definite demands that A + /3e m i n > 0. Consequently we 
have 

A > -/3e min ■ (10) 

Taking an AFM Jx, we consider two special cases. First, 
for an AFM J 3 and all other J n >i = 0, we have e m i„ = 
— I6J3. Second, for a FM J 2 and all other J n >i = 0, the 
position of the minimum in reciprocal space is given by 
Chern et al— in their Eq. (4), with e min ~ — 8.57| J 2 | for 
I J 2 | « Ji- 

III. CONSEQUENCES OF AFM J 3 OR FM J 2 

A. Paramagnetic phase 

Within the paramagnetic phase, AFM J3 or FM J 2 
act to suppress pinch points, and to promote cluster-like 
correlations. These effects are shown in Fig. [31 which 
should be compared, for example, with Fig. 3 of Ref. Q. 
The behaviour illustrated here is one of our main results. 
Its origin is discussed in this section. 

1. Consequences for pinch point visibility 

The pinch points are the signature of algebraic cor- 
relations in a Coulomb phase. There are three ways in 
which they may become less visible, (i) Further neigh- 
bour interactions or finite temperature endow the pinch 
points with a width set by the inverse of the longitudinal 
correlation length, (ii) The absolute amplitude of the 
pinch points is set by the inverse of the flux field stiff- 
ness, (iii) The visibility of the pinch points depends also 
on the relative weight at the pinch points compared to 
that at the maximum in scattering. We address these 
mechanisms quantitatively below, (i) For AFM J3, the 
longitudinal correlation length is ay/ f3(Jx — %J^)/\. The 
lower bound on A puts an upper bound on the correlation 
length of a v / {J 1 /16J 3 - 1/2). Even with J 3 ~ J x /150, 
the correlation length never exceeds 3 site spacings and 
the pinch points are thus never sharp. For FM J 2 , a 
similar argument gives an upper bound on the correla- 
tion length of approximately ay/ (2 + J\/\ J 2 |)/8.6. The 
effect on the correlation length is shown in Fig. Ufa) for 
J 2 = —0.025 J\. (ii) Since the absolute amplitude of the 
pinch points is set by 1/A, the lower bound on A puts an 
upper bound on the absolute intensity of pinch point cor- 
relations. For AFM J3, recall that A > 16/3 J3, whereas 



for FM J 2 , we have A > 8.6/3J 2 , for | J 2 | < Jx- The pinch 
point intensity is shown in Fig. Hl^b) for J 2 = — 0.025 J\. 
(iii) Moreover, as spectral weight shifts away from the 
pinch points to other wavevectors, the relative amplitude 
given by the ratio of pinch point correlations to the max- 
imum in the correlation function decreases even more 
rapidly. These three mechanisms indicate how pinch 
points may both be suppressed in amplitude and broad- 
ened in width by further neighbour interactions. 



2. Hexagonal cluster-like correlations 

We have seen how further neighbour terms suppress 
spectral weight at the pinch points. We turn now to why 
antiferromagnetically correlated loops of spins such as 
hexagonal clusters may appear in a model with further 
neighbour terms. The mechanism is straightforward to 
explain in the long wavelength description. For simplic- 
ity, we consider first the limit (3Jx — > 00 (at finite /3J 2 ). 
In this case, pinch points remain infinitely sharp since the 
local constraint of zero divergence is enforced rigorously, 
but their amplitude is still suppressed since the flow of 
the stiffness with temperature is unaffected. The further 
neighbour terms couple to transverse gradients of the flux 
fields and the effect of AFM J3 or FM J 2 is to promote 
transverse gradients in the flux field. In the continuum 
picture, maximising transverse gradients requires having 
as many small loops of flux as possible; when understood 
in terms of the lattice geometry, the smallest loops of 
flux available are the hexagonal clusters. Considered mi- 
croscopically, it is also clear from Fig. [T] that an AFM 
Jsb is likely to induce antiferromagnetically correlated 
hexagons of spins. 

Following Ref. 0, we define Hq as the temperature- 
dependent position of the maximum in scattering along 
the (hh2) direction (see dashed lines in Fig. [5]). For FM 
J 2 , at T=0, Hq coincides with h* given by Chern et al. 19 
in their Eq. (4), the position of the minimum eigenvalue 
of the exchange matrix. The evolution of Hq with tem- 
perature extracted from the SCGA is shown in Fig.@|c). 
Qualitatively similar behaviour for Hq is found experi- 
mentally in CdFe 2 04 (see Fig. 4 in Ref. d). In contrast, 
in the NN model (at least classically) the scattering max- 
ima are always found at the pinch points, i.e Hq — 0, so 
the NN model does not provide an explanation for why 
the maximum in scattering is elsewhere. 

In Fig. |6j we show the structure factor on the same 
planes of reciprocal space as in Fig. [3J but for the case of 
a FM J 2 rather than an AFM J3. The qualitative effect 
is similar, but the characteristic scattering shapes look 
less like the hexagonal form factor used to explain the 
scattering in many chromites. 
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FIG. 3. (color online) Structure factor S(q) at T = O.lJi. (a) and (b): NN interactions only; (c) and (d): J3 = 0.025Ji. (a) 
and (c): (hOl) plane; (b) and (d): (hhl) plane. Contour lines are at the levels marked on the color scale. Compared to a nearest 
neighbour model, the effect of AFM J3 is, first, to suppress and broaden the pinch points, and second, to shift spectral weight 
so that behaviour resembles more the hexagonal cluster scattering. MC simulations yield plots that are indistinguishable by 
eye. 




FIG. 4. Temperature dependence of quantities characterising 
the scattering, for J2 = — 0.025 Ji. (a) Coulomb phase corre- 
lation length £ in units of site spacing for NN model (broken 
line) and for model with non-zero J2 (solid line), (b) Ampli- 
tude of correlation function at pinch-point, for both models, 
(c) Position Hq of maximum in scattering, as described in 
text. Dashed line indicates temperature T = J2: the effect 
of non-zero J2 extends to T > J2. Curves are from SCGA. 
Points are from simulations on 2048 spins. 



B. Ordering transitions 

For J 3 <C Ji, we observe a first order phase transition 
at T c ~ 2.2 J3 in our MC simulations (eg. see upper inset 
in Fig. [5]). We have not attempted to characterize the or- 
dered state since our focus is on the strongly correlated 
paramagnet in the window T c < T < J\. A phase tran- 
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(hh2) (r.l.u) 

FIG. 5. Structure factor along (hh2) for J2 = — 0.025 Ji. Lines 
are from SCGA, points from simulations on 16384 spins. The 
position Hq at each temperature is indicated by dashed line. 



sition is not predicted by mean field calculations^, nor 
within the theory of Sec. Ill Al since the interaction min- 
ima has a 1-dimensional degeneracy in reciprocal space. 

In the phase diagram presented by Chern et al.— , fer- 
romagnetic J2 leads to incommensurate order at low tem- 
peratures consistent with mean field calculations^, but 
also an intermediate collinear regime, stabilised by ther- 
mal fluctuations, with first order transitions out of the 
paramagnetic phase. The transition temperature is ap- 
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FIG. 6. (color online) Structure factor at T = O.lJi with 
J 2 = -0.025 Ji. Panels as in Fig. [3] 



proximately at T c ~ J2, but see Ref. [l_9| for details. 



IV. CONSEQUENCES OF AFM J 2 OR FM J 3 

Although the main focus of this paper is further neigh- 
bour terms with their sign chosen to suppress pinch point 
scattering, we now also consider further neighbour terms 
with the opposite sign. These enhance pinch point scat- 
tering. 



A. Paramagnetic phase 

It is clear from Eq. ([6]) that in the flux description an 
antiferromagnetic J2 or FM J3 penalizes transverse gra- 
dients of the flux fields. Qualitatively then, one expects 
spectral weight to accumulate at small q as short wave- 
length fluctuations are suppressed, developing into Bragg 
peaks as the transverse flux correlation length diverges at 
a phase transition. 

As seen from the numerical solution of Eq. © shown 
in Fig. [5J FM J 3 causes the stiffness A to decrease with 
cooling, and the same behaviour is found with AFM 
Ji- The strength of flux fluctuations is given by the 
inverse stiffness as in Eq. ([7]) and consequently this in- 
creases. In Fig. [7] we show spin correlations in the para- 
magnetic phase at the same value of (3Ji as in Fig. [6] 
with antiferromagnetic J2, pinch points are intensified 
and sharper than for the nearest neighbour model as ex- 
pected. Youngblood and Axe— considered the Maxwell 
action with terms penalising transverse flux gradients 
and obtained a similar effect. 
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FIG. 7. (color online) Structure factor S(q) at T = O.lJi 
and J2 = 0.025Ji in (a) the (hOl) plane and (b) the (hhl) 
plane. Contour lines are at the levels given in the color scale. 
Positive J2 increases pinch point scattering, and eventually 
leads to Neel order. The shapes around the pinch points are 
almost identical to those discussed by Youngblood and Axe—. 



B. Ordering transitions 

With an AFM J 2 < Ji, the full model, Eq. (P, has 
a first order transition to a q = Neel-ordered phase 
at T ~ 3.2J 2 19 . Within the framework of Sec. UTaI the 
transition is second order and occurs at a lower temper- 
ature (T c ~ 2.44J 2 for J1/J2 = 100): the one at which 
the minimum eigenvalue of the quadratic form in Eq. @ 
vanishes. For FM J3 -C Ji, our simulations reveal a first 
order transition at T ~ S | J3 1 (see lower inset of Fig. [2]). 
In both cases, the transition is to a fully flux-polarised 
state. Similar transitions occur in other Coulomb phase 
systems: for pyrochlore Heisenberg spins, this is q = 
Neel order, for interacting dimer models, a transition to 
a uniform tilted phase, and in the context of polarization 
fluctuations in paraelectrics^, a ferroelectric transition 
out of the paraelectric phase. 

The transverse flux correlation length £j_ , which is zero 
for a pure nearest neighbour model, diverges as A — > 0. 
Flux fluctuations are correlated over larger regions and 
the divergence of £_l marks the transition into an ordered 
phase. 



V. GENERAL COMBINATIONS OF FURTHER 
NEIGHBOUR INTERACTION 

In Sees. Mil and IIVI we took both third neighbour 
couplings as equal, but in general they need not be 
and here we consider J^ a and J$b as distinct. Within 
the long wavelength description, nothing qualitatively 
changes since further neighbour terms can only act to add 
gradient terms while renormalising the stiffness. How- 
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FIG. 8. S(q) in the (hhl) plane for h € [0,3/2], I € [1,3] 
at T = O.lJi, J2 = and with Js a , J$b interactions treated 
independently, showing the variety of short range correlations, 
and the effects on the pinch point at (002). 



co=0.75 J /, 

1 f/ 




FIG. 9. Points: ^(q,^) along (hh2) at T = O.lJi, J 3 = 
0.025 Ji from simulations. Solid line: S(q, oS) from stochastic 
model of Ref. [2?] with further neighbour terms. Broken line: 
S'(q) rescaled for comparison. Pinch points are more sup- 
pressed at small-o;, and conversely less suppressed at high-o;. 



ever, the detailed nature of the short range paramagnetic 
correlations is not generic and depends sensitively on the 
ratios of the different further neighbour couplings. Since 
this balancing act is separate for each candidate material 
which our ideas hope to describe, we provide a survey for 
several different combinations. In Fig.[8j we have plotted 
paramagnetic correlations in the reciprocal space patch 
described in the caption for a variety of different J^ a and 
J 3 b at T = O.lJi. By the arguments of Ref. [T^, J 2 acts 
as —Jza for large enough J\, so the horizontal axis can 
be considered approximately to probe J^, a — Ji- 



VI. DYNAMICS AND INELASTIC 
SCATTERING 

We have so far considered only equal time correlations 
but further neighbour (FN) terms will also affect dynam- 
ics. There is an intrinsic precessional dynamics associ- 
ated with the Heisenberg spins in Eq. ([T]). Nevertheless, 
we have shown elsewhere^ for the NN model that the 
low temperature dynamics is dominated by relaxational 
modes and is very well captured by a purely relaxational 
stochastic model based on the SCGA of Sec. Ill Al We 
find this statement about the low temperature paramag- 
net remains true even with FN terms. As evidence for 
this, we present in Fig. |H] the results of molecular dj 
namics simulations using the same methods as in Ref. 
for the dynamic structure factor, plotted at 3 different 
frequencies co. On the same plot, we show the predic- 
tion obtained from including FN terms in the stochastic 
model. Compared to the static result, the pinch point 
scattering at (002) is more suppressed at small to, and 



f.S 



conversely, less suppressed at higher to. The implications 
for quasi-elastic scattering experiments at to ~ are that 
we expect the effects of FN terms to be more conspicuous 
(and conversely at larger to, to become less conspicuous) 
when compared to the static (w-integrated) scattering 
which is the main focus of this paper. 

The enhanced suppression can be understood by con- 
sidering the stochastic model in Ref. l27|. Within that 
model the dynamic structure factor is understood as the 
sum of contributions from the four bands of the interac- 
tion matrix, given by 



4 

Sfaw) = ^^(q)/ (U (q, 

,, — 1 



to 



(11) 



where the <7/t(q) are form factors, independent of tem- 
perature, and where the thermal weight and dynamics 
for each band are contained in the functions / M (q, to). 
The quadratic form in Eq. @ has eigenvalues a M (q) , la- 
beled by band index fj,, which determine relaxation rates 
r^(q) oc cv M (q) and these enter the dynamic structure 
factor through the relation 



Ufa,") 



^(q) 



2iyq) 

2 + iVq) 



fl2) 



By integrating over lu we obtain the static structure fac- 
tor, which has the form of Eq. (lTTj) . but with a~ 1 (q) in 
place of / M (q, to). 

With only nearest neighbour interactions, the two 
lowest bands (which are dominant at low temperature) 
are dispersionless and degenerate, which has two conse- 
quences: the q-dependence arises entirely from the form 
factors in Eq. since the relevant are independent 
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of q, and the w-dependence approximately factorises out. 
With FN interactions, the lowest bands acquire disper- 
sion and the relevant are no longer independent of 
q. If their ratio in the static case for band \i at qi and 
q2 is r, then in the dynamic case with ui = it is r 2 . 
This is why the suppression effects we have discussed in 
this paper, which correspond to r 1, arc enhanced in 
quasi-elastic scattering. Further discussion of dynamical 
correlations is given in Ref. [2l|. 



VII. ESTIMATES OF FURTHER NEIGHBOUR 
EXCHANGE 

Evidence from ab initio calculations suggests third 
neighbour couplings are important in a number of rele- 
vant compounds. In a study of ACT2X4 spinels including 
zinc and cadmium chromite, Yaresko™ finds that AFM 
coupling between third neighbours is important in all the 
compounds considered, and that the ratio Jsb/Jza — 0.5 
for Zn and Cd. (Our J^ a and J^b are labeled as J& 
and Js'i in Ref. |28[ ). A separate ab initio study of 
CdCr 2 04 finds J3/J1 ~ 0.3 where both third neighbour 
couplings are taken with the same strength 2 ^. In both 
ZnFe2C>4 and CdFe204, Cheng— also reports that third 
neighbour interactions are estimated to be much stronger 
than second-neighbour interactions and on the same or- 
der in magnitude as first-neighbour interactions. This 
result is in broad agreement with a study of third neigh- 
bour interactions in ZnFe2C>4 ( Ref. f20h . 

Experimental estimates of nearest neighbour and 
next nearest neighbour exchange from susceptibility 
measurements 3 - , find for CdCr 2 04 a FM J 2 ~ — 0.3 Ji 
but for ZnCr 2 4 an AFM J 2 ~ 0.15Ji. From our re- 
sults in Sec lIVl the latter would imply scattering quite 
different from that observed; we have not examined the 
reliability of the model used in Ref. in the presence of 
third neighbour terms, but believe that ZnCr204 is un- 
likely to be described by a Heisenberg model with AFM 
J2 without other further neighbour terms included. 



VIII. DISCUSSION 

The appearance of antiferromagnetically correlated 
hexagonal clusters in the spinels is certainly puzzling 
in the context of a pure nearest neighbour model, since 
properties of the NN model have been reliably calculated 
and give different behaviour. In Sec. IVII1 we discussed 
evidence for further neighbour exchange and found that 
AFM third neighbour exchange is important in many of 
the compounds we have been discussing. We should also 
address other possibilities. 

It has been suggested that spin clusters are induced by 
magnetoelastic coupling 3 . However, we note that among 
the compounds where spin clusters have been reported, 
CdCr 2 04 undergoes a c-axis elongation^, whereas both 
ZnCr204 and MgCr204 undergo a c-axis contraction, 



and CdFe204 has no structural transition at all; it is not 
clear how magnetoelastic coupling which leads to differ- 
ent structural transitions should explain the cluster-like 
scattering common to all of these materials. Moreover, 
scattering in Y(Sc)Mn2 is broadly similar to the spinels 6 , 
but we expect the nature of spin lattice coupling to be 
different, and this throws further doubt over the idea that 
magnetoelastic coupling may be universally responsible 
for cluster-like scattering. 

We should also discuss the temperature dependence of 
our results. The images in Figs. [3] and [5] are for a sin- 
gle temperature, where the hexagonal cluster-like scat- 
tering is most visible; at lower temperatures however, 
the further neighbour terms induce magnetic phase tran- 
sitions as discussed in the text. If indeed further neigh- 
bour terms are required to explain the hexagonal cluster 
scattering, we may also expect to observe the associated 
magnetic phase transitions. In the chromites, we can 
avoid the issue by assuming the magnetic phase tran- 
sition is preempted by the structural transitions driven 
by the independent mechanism of spin-lattice coupling; 
this is not entirely satisfactory as an explanation, since 
CdFe 2 04 shows neither a structural transition, nor devel- 
opment of any Bragg peaks at any temperature. It is pos- 
sible that quantum fluctuations or disorder can stabilise 
the paramagnetic phased, but it is hard to be conclusive 
about a range of materials at once, and the role of fur- 
ther neighbour exchange should be assessed individually 
for different materials. 

We have presented further neighbour exchange as the 
simplest extra interaction which gives hexagonal cluster 
scattering, but further neighbour exchange is not the 
only perturbation one could imagine. Others include 
Dzyaloshinsky-Moriya interaction s 32 ' 33 and biquadratic 
exchange, which can appear microscopically or as an ef- 
fective term taking into account quantum fluctuations^. 
In future work it would be interesting to consider the ef- 
fects of these other perturbations on the paramagnetic 
phase. The present results are sufficient to show that 
further neighbour exchange of a realistic strength can 
dramatically alter the diffuse scattering and should not 
be ignored in comparing experimental results to simple 
theoretical models. 

Clusters should be understood, then, as short range 
order in a strongly correlated state. A similar view is 
expressed by Yavors'kii et al— when considering cluster- 
like scattering in spin ices; their analysis including further 
neighbour terms suggests that cluster-like scattering is 
the property of a strongly correlated liquid state which 
is sensitive to weak perturbations rather than due to the 
emergence of 'real' clusters. 

Our results potentially explain why a clear signature 
of algebraic correlations is missing in candidates for py- 
rochlore Heisenberg magnets. Small further neighbour 
interactions are sufficient to wash out the pinch points. 
Our results may have implications for other pyrochlore 
models, but in Ising models where the ice-rules are expo- 
nentially enforced, we do not expect the dramatic broad- 
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ening we observe in the Heisenberg model. Nevertheless, 
the suppression of pinch point scattering amplitude is ex- 
pected, and we note that in Ref. d such a suppression is 
reported in the spin ice Ho 2 Ti 2 07. 

In summary, we have studied the effect of further 
neighbour interactions on the low temperature para- 
magnetic phase of the frustrated Heisenberg antifcrro- 
magnet on the pyrochlore lattice. Further neighbour 
terms induce transitions to ordered phases, but they also 
have a striking effect on paramagnetic correlations. In 
the description of the low temperature paramagnet as 
a Coulomb phase, further neighbour terms cause the 
Coulomb phase coupling constant, or stiffness, to flow. 
Pinch points in diffuse scattering have their amplitude 
and width controlled by the stiffness, and even very weak 
further neighbour terms can cause pinch points to be sup- 
pressed in amplitude and broadened. With FM J 2 or 
AFM J3, paramagnetic scattering is altered to resemble 
more the hexagonal cluster scattering that is often ob- 
served in experimental systems, notably spinels. Further 
neighbour terms then provide a mechanism for the pre- 
viously unexplained cluster-like scattering in frustrated 
spinels. 

ACKNOWLEDGMENTS 

PC acknowledges helpful discussions with M. Gingras 
and S. Bramwell. This work was supported in part by 
EPSRC Grant Number EP/D050952/1. 



where r runs over N fee lattice sites so that there are 4N 
spins in total. 



Appendix B: Interaction matrices 

As described in the main text, with A ™ denoting the 
adjacency matrix that is 1 if sites i and j are nth neigh- 
bours and zero otherwise, the interaction matrices are 
chosen to be v£ } = A$ + 2% V% ] = Af> + 4%; 

Vff a) = and V^ b) = A<? b) -6dy. This choice 

of diagonal entries is made so that when expressed in re- 
ciprocal space according to Eq. (|A3|) and using the basis 
change in Eq. ((5]), the 3x3 subblock which acts among 
the flux components vanishes at q = 0. This ensures that 
A can be interpreted as the stiffness of the flux fluctua- 
tions. 

The matrix elements for the interactions can be found 
in Ref. fvK although in a different basis so we provide 
them here using our conventions. We first define the fee 
displacement vectors — 2(c M — c„) where the c M are 
as in Appendix. [21 After Fourier transform the matrix 
elements of the adjacency matrices are given by 



^ (1) (q)i2 =2 




and 



Appendix A: Lattice definitions 

We take throughout a conventional cubic unit cell of 
side 1, and use a to refer to the pyrochlore nearest neigh- 
bour distance, so that a 2 — 1/8. The midpoints of 
th e tetr ahedra form a diamond lattice with site spacing 

The pyrochlore lattice is a face centered cubic (fee) lat- 
tice decorated with a tetrahedron at each fee site. Using 
the following choice of fee basis vectors 



ai 



; {e y + e z ), a 2 = -(( 



■ z +e x ) , a 3 = -(e x + e y ) , 
(Al) 

then (up to an arbitrary displacement) the corners of the 
tetrahedra appear at positions 



In — — f^A 

I 2 '2 ' 2 J 



(A2) 



where /i labels the four sublattices and runs from 1—4, 
as in Fig. [TJ The Fourier transforms we use are 



^« = -7= E ^(qK 

^ qGB.Z 



q.(r+c^) 



'^) = 7=E^(r)e- i;q - (r+c 



(A3) 



,(2)/ x o (tl3 +t 23 )\ . „ /q • (tl4 +t 24 ) S 

A> A} (q)i2 = 2 cos I I +2 cos I 



with diagonal elements zero, and the rest by permutation. 

Furthermore, third neighbour couplings are diagonal 
giving: 

-4 (3a) (q)n =2cos(q-t 12 )+2cos(q-t 13 )+2cos(q-t 14 ) 
and 

= 2 cos (q • t 2 3)+2 cos (q • t 24 )+2 cos (q • t 34 ) . 

To derive the long wavelength theory in the flux sector, 
we implement a change of basis as described in the main 
text using the orthogonal matrix 



/l 1 1 1 




and then take the small q limit. This procedure gives the 
following interactions 
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and 



p V (2) P T^ 



and 



/8 



PV^P 1 



a 2 g 2 



-a"q x q z 

-,2 



-a 2 q y q z 

a 2 q 2 x 

a 2 q x q y 



\-a z q x q y a 2 q x q z 



-a 2 q x q z ~a 2 q x q y \^ 

a 2 q x q y a 2 q x q z 

a 2 q 2 a 2 q y q z 

i ii 
a q y q z a z q; 
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2a 2 q y q z 
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\ 2a 2 q x q v 



2a 2 q y q z 
-2a 2 q 2 x 
-2a 2 q x q y 



2a 2 q x q z 2a 2 q x q y \ /O 



-2a 2 q x q y -2a 2 q x q 
— ,,,, ~2a 2 g 2 -2a 2 q y q 
-2a 2 q x q z -2a 



,2 
2 

-2a 2 q 2 ) 



-4a 2 



qi + q 2 z 












PV (3a) P T ~ 4a 2 
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~q x q z 

\-q x q y 



-q y qz 
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q x q y 
-q x q z 



-q x q z 
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-q 2 y 
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-q x q y \ 

-q x q. 



-qyqz 
-q 2 z 



/o 



4a 2 



qi 









+ g: 
g 2 









3* 






i 



2 

%) 



\ 



PV ( 3b )p T ~4a 2 



/-g 2 
q y qz 
q x q z 

\q x q y 



q y q z 
ql 

q x q y 
q x q z 



q x q z 
q x q y 

q 2 
hi 

q y q z 



q x q y \ 

q x q z 
q y q z 

ql 



4a< 



( 2 

2g 2 




q 2 












q 2 y 



Hi) 



In each case we have split the interaction into a part 
which acts like J\ within the flux field subspace, and 
another part. Chern et al— argue that J 2 and — J 3a 
have the same effect if Ji — > oo. This fact is evident 
in the small q expansions above. If the two symmetry 
inequivalent third neighbour couplings are taken to have 
the same strength, as done by Reimers et al.—, then the 
rotated third neighbour interaction matrix is particularly 
simple: 

P[V< 3a > + V( 3b >]P T ~ -8a 2 diag(g 2 ) . (Bl) 

Integrating out M(q) in the long wavelength limit leads 
to couplings at C(g 4 ) between the components -Bj(q), 
which we have ignored in the analysis at small q in 

Sec. urn 

Appendix C: Locations of pinch points in structure 
factor 

The experimental probe in scattering experiments is 
the structure factor 5(q) = 4(|M(q)| 2 ). However the 
long- wavelength physics is in the field B near q = 0. In 
this appendix we determine the places in reciprocal space 
where signatures of the behaviour of B are visible. 

For K a reciprocal lattice vector we have 

5 A1 (q + K) = e 4K - c -5 A1 (q). (CI) 
Taking as the reciprocal basis to slj , then a reciprocal 



lattice vector K = nibi + n 2 hi + n^hi is 

K = 2n[(n 2 + n 3 — ni)x 

+ (n 3 + n x - n 2 )y + {n x + n 2 - n 3 )z] (C2) 

Bearing in mind the definitions of c^, the phase factors 
for the different sublattices are 1, (-1)™ 1 , (-1)" 2 , (-1)™ 3 . 
It follows that 

M(q + K) = i[«oM(q) + vB(q)] (C3) 

where 

t; = l + (-l)" 1 + (-l)" 2 + (-l) 113 (C4) 

and 

/l + (-1)" 1 - (-1)" 2 - (-1)" 3 \ 
v = 1 (-1)" 1 + (-1)™ 2 - (-1)" 3 . (C5) 

yi - (-i) ni - (-i)" 2 + (-i)™ 3 / 

If all rii are even, the structure factor in the vicinity of 
the reciprocal lattice point probes only the total magne- 
tization. If only one rii is even, it probes v • B, and in the 
two other cases, there are contributions from both M(q) 
and v • B(q). 

Appendix D: Other results 

We also provide here the structure factor for the near- 
est neighbour model. In the notation of Isakov et al.—, 
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where c a b 

Q c xy 



,, 2 
■'.V 



and 



r» 2 



a 6 



cos ( qa 4 qb ) with defining s 2 , = sin 2 (^-) and 
r - 3 and further structure factor is 

I 



■(ab) 



Cab 



c-y, the static 



S(q), 



(2 - C(^) - c {yz) - c {zy) )t\ + \{t\) 2 - 4 [c^s 2 + C( yz) s 2 x + c (za ).' 
3-Q + 2tA+ i(iA) 2 



(Dl) 



where t = T j J\ and A is the solution of Eq. ([3]). At due to the choice of spin length, Eq. (|D1|) reduces in the 
T <C Ji, the zero temperature solution A = 3/2 is a good limit t — > to the zero temperature form provided by 
approximation (see Fig. [2]). Up to an overall rescaling Isakov et al.—. 
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